Theory of warm ionized gases: equation of state and kinetic Schottky anomaly 
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Based on accurate Lennard- Jones type interaction potentials, we derive a closed set of state 
equations for the description of warm atomic gases in the presence of ionization processes. The 
specific heat is predicted to exhibit peaks in correspondence to single and multiple ionizations. Such 
kinetic analogue in atomic gases of the Schottky anomaly in solids is enhanced at intermediate and 
low atomic densities. The case of adiabatic compression of noble gases is analyzed in detail and the 
implications on sonoluminescence are discussed. 
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The recent progresses in stellar astrophysics 
experiments on inertial confinement in plasmas 
and in plasma diagnostics have triggered renewed inter- 
est in the determination of accurate state equations for 
mixed systems of atoms and ions in ionized gases 043- 
At the same time, much recent effort has been devoted 
to investigating the role of plasmas of ionized atomic 
gases in sonoluminescence (SL) 0|-18|. Sonolumines- 



cence consists in the emission of short flashes of light 
from collapsing bubbles in a liquid excited by ultrasonic 
waves. The light-emitting bubbles contain mainly noble 
gas atoms plus small quantities of water vapor and in the 
last part of the collapse they can reach densities compa- 
rable to those of solids. According to recent results [18| 
the conditions inside the bubbles may be very extreme, 
the temperature may far exceed those estimated from 
the emitted light, and the plasma electron density can 
reach the same order of magnitude as that of densities 
created in laser-driven fusion experiments. A detailed 
understanding of the noble-gas thermodynamics is thus 
crucial in order to explain the experimental results and 
uncover the basic mechanisms of sonoluminescence. 

In the present work we analyze the thermodynamics 
of atomic gases for different regimes of densities and 
temperatures characteristic of warm plasmas, and con- 
sider applications to the case of noble gases by studying 
regimes that can be typical of the last stage of collaps- 
ing bubbles in sonoluminescence, i.e. the stage of light 
emission. We derive a set of state equations based on ac- 
curate microscopic interaction potentials of the Lennard- 
Jones type and we predict the presence of anomalous 
peaks in the specific heats. These peaks are explained 
in terms of a kinetic extension of the Schottky anomaly 
observed in solids with a limited number of energy lev- 
els at very low temperature fl9| . The peaks appear in 
correspondence to ionization processes as the tempera- 
ture is increased. This implies that the different degrees 
of ionization play in plasma physics a role similar to the 
one played by the finite number of energy levels in the 
solid state case, the fundamental difference lying in the 
fact that the predicted nontrivial contribution to the spe- 



cific heat in ionized gases is provided by the unbounded 
nature of the atomic kinetic energy. The existence of ki- 
netic Schottky anomalies bears important consequences 
on the collective heating of atomic gases. In particular, 
collective heating is strongly suppressed when the spe- 
cific heats undergo a sharp increase in correspondence of 
atomic ionizations. We will discuss in detail the case of 
noble gases, which is especially relevant for sonolumines- 
cence and various interstellar phenomena. Considering 
the late stage of adiabatic compression, we determine 
self-consistently the number density of plasma electrons 
and find good agreement with recently reported experi- 
mental results [18j |. 

Equation of state - By taking into account the kinetic, 
interaction, and ionization contributions, the energy den- 
sity of the atomic ensemble can be expressed as 
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where the index i classifies the particles present in the 
systems (electrons, atoms, and ions of different charge 
multiplicity) by their charge expressed in units of |e|. 
Explicitly, i — — 1 for free electrons, i — for neutral 
atom and j g [1, Na] for ions of different charge, rang- 
ing up to the maximum allowed value corresponding to 
the atomic number Na- The number density of the i— th 
particle species is denoted by rij, while Ei is the energy 
needed to ionize an atom i times, with E—\ = Eq = 0. Fi- 
nally, hi is the potential energy density of the i-th particle 
species due to the interaction with the other particles of 
the system. Taking the virial expansion up to second or- 
der (higher orders may be considered when necessary), 
hi takes the form 
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where Uij(r) is the interaction energy between the i-th 
and the j-th particles at distance r. 

The atom and ion densities n, are not fixed a priori, 
but are determined by the underlying quantum dynamics 
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through the self-consistent coupling of Eq. ([5J with the 
Saha equations [20j : 
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Here m e is the electron mass and gi is the ground-state 
degeneracy of an atom ionized z-th times. Eq. ([3]) is a set 
of Na equations in Na + 2 variables represented by the 
densities rij. The set can be closed taking into account 
that the electronic density is related to the ion density 
by 7i_i = ™! an d that no is related to the total 

atomic density n by no = n — X)<=a n i- 

We will now provide a closed expression for the pres- 
sure P. In agreement with Eq. (TTJ), P needs to be written 
as the sum of a kinetic term and a term associated to the 
interaction of every particle with the surrounding ones, 



N A 

P = K b T 



2tt 



N A 

E 



(4) 



where the interaction term /ij, considering the virial ex- 
pansion up to second order, reads 
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By means of Eq. (fTJ) wc can compute the specific heat at 
constant volume, c v — (§f^) y . Eq. <(4]) can be used to 
obtain the specific heat at constant pressure by using the 
'ar^ -1 (dP\ 2 



relation 



-T 



§£) T {%) v . Eqs.©-© iden- 
tify a complete set of state equations for a gas in which 
the degree of ionization is determined sclf-consistently 
from thermal equilibrium. 

Non interacting systems: warm plasmas at low density 
- We consider first the situation of ultra-low densities, 
such that the interaction effects can be neglected, i.e. we 
can put (Uij = 0, Vi,j). For this regime, the main 
results are illustrated in Figs. [T] and [2] for a system of He 
atoms. We investigate the behavior of the pressure P, 
the specific heats c v and c p , and the density of free elec- 
trons n_i as functions of the temperature T at different 
values of the total density. The top panel of both figures 
refers to the instance of n = I0 10 cto~ 3 , the typical den- 
sity in the active galactic nucleus and about ten orders 
of magnitude less than the corresponding value at room 
temperature on earth's surface. From Fig. [1] wc observe 
that both specific heats exhibit two peaks in exact cor- 
respondence, respectively, with the processes of first and 
second ionization, as identified by the associated jumps 
in the free-electron density rt_i. In Fig. [2] the pressure 
shows a step-wise behavior corresponding to the increase 
in the number of free electrons released in the ionization 
process. Going to regimes of much higher densities, as 
reported in the central and bottom panels of Figs. [1] and 
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FIG. 1: The molar specific heats at constant volume c v (Red 
solid line), and at constant pressure c p (black dot-dashed line) 
as functions of the temperature T, compared with the electron 
density n_i (blue dashed line) for an ensemble of He atoms 
at three different densities n. Top panel: n = 10 10 cm _ ' ! ; 
central panel: n = 10 19 cm -3 ; bottom panel: n = 10 22 cm -3 . 
The specific heats are expressed in units of the universal gas 
constant R. 



[2 we observe the progressive smoothing out of the peaks 
in the specific heats, as well as the broadening of the step- 
wise behaviour of the density n_i and of the pressure P. 
This is due to the fact that as the density increases, the 
respective ranges of temperature in which the processes of 
first and second ionization occur start to merge. The on- 
set of the peaks occurs as the temperature increases and 
the corresponding thermal energy approaches the differ- 
ence between the atomic energy levels. In this situation 
the different ionization states become significantly pop- 
ulated and hence there is an enormous increase in the 
entropy of the system due to the uncertainty on the de- 
gree of ionization. This implies the presence of peaks 
in the profiles of the specific heats. At variance with the 
Schottky anomaly occurring in solid state physics for sys- 
tems with a finite number of energy levels (bound states), 
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FIG. 2: The pressure P (red line) as function of the tempera- 
ture T, compared with the electron density n_i (blue dashed 
line) and with the pressure obtained from the ideal gas law 
neglecting ionization processes (black dot-dashed line), for an 
ensemble of He atoms at the same three values of the density 
reported in Fig. \T\ 



in the case of non interacting plasmas the rate of change 
of the entropy is due to the kinetic energy of the ionized 
atoms, resulting in multiple peaks of the specific heats, 
one for each successive degree of ionization. Therefore 
this a kind of kinetic Shottky anomaly, in principle with 
an unbounded number of peaks in the specific heats. 

Interacting systems: warm dense plasma - We now 
include the interactions between different particles in 
our investigation by introducing two-point interatomic 
potentials f/jj which can describe with reasonable ac- 
curacy the microscopic physics of plasmas in the in- 
termediate regimes of temperature and density, i.e. 
warm dense plasmas. Here and in the following 
we will model the interaction between two neutral 
atoms of a noble gas by Lcnnard- Jones type poten- 
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in which e stands for 



the depth of the potential well and a is the finite dis- 



tance at which the inter-particle potential vanishes. The 
Lennard- Jones potential that we propose to adopt is 
well justified as it amounts to the sum of a short-range 
Pauli repulsion due to the partial overlap of the elec- 
tronic clouds, and a long-range dipole-dipole attraction 
between the two neutral atoms. We need also to in- 
clude the interaction energy between a neutral atom and 
an ion, which can be described by adding to i/o,o the 
chargc-dipole attractions Uo.i(r) = Uo,o(r) — i 2 aEf(r), 
where a is the polarizability of the neutral atom and 
Ei(r) is the electric field associated to the electrostatic 
potential Wi(r) that solves the corresponding Debye- 
Hiickel equation [20j . Following the same idea, we de- 
scribe ion-ion interactions by further adding the chargc- 
dipole interactions and the electrostatic ion-ion repulsion: 
Ui,j(r) = Uo,o(r) - (* 2 + j 2 )aEf (r) + ijWi(r). 

The interactions involving one or more free electrons 
can be described by considering any free electron as a 
cloud similar to the one surrounding an Hydrogen atom, 
but with a radius equal to the thermal De Broglie wave- 
length A and a vanishing polarizability. With these as- 
sumptions, free electrons can be treated as ions and their 
interactions similarly to the ones between positive, mas- 
sive ions. We have only to replace, in the Lennard Jones 
term, a with ^-j^ and e with y/e£_H (being eh the depth 
of the potential well of the Lennard- Jones potentials be- 
tween two Hydrogen atoms) in the case of an interaction 
between a free electron and an ion, and to replace a with 
A and e with ejj in the E/_i._i(r) term. 

In order to compare the cases with and without inter- 
actions, we consider here and in the following a system of 
Xenon atoms at high density. The choice of Xenon is mo- 
tivated by the fact that this element is characterized by a 
lower ionization energy than that of Helium [2l[. Hence, 
one may reasonably expect that, especially at high den- 
sity, the corrections due to interactions will be sizeable 
and clearly detectable. In our further treatment we will 
neglect the dependence of a and e on the degree of ion- 
ization. This approximation does not affect results in the 
case in which one considers only the first few ionizations. 
Therefore, in the analysis summarized in Fig. [3] we take 
into account only the processes of first and second ion- 
ization in Xenon at densities of the order n = 10 22 cm -3 . 
Fig. [3] shows that particle interactions do not modify sig- 
nificantly the kinetic Shottky anomaly in c v , even at high 
densities. Indeed, the fields hi in Eq. fl2J) are about one 
order of magnitude lower than the ionization energy and 
therefore their effect on the solution of the Saha Equation 
Eq. ([3]) is negligible. On the other hand, the behaviors 
of P and of c p are strongly influenced by the presence 
of interactions, as the latter introduce relevant Newto- 
nian effects in the conservative forces and in the ensuing 
associated pressures, whose net effect is the strong sup- 
pression of the Schottky-likc type of behavior. 

Kinetic Schottky anomaly and adiabatic compression 
- We now discuss the influence of the kinetic Schot- 
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FIG. 3: Behavior of the specific heats as functions of the tem- 
perature, with and without interactions, for a molar weight 
of Xenon. Main plot: c v for the non interacting case (dot- 
ted black line), c v for the interacting case (solid red line), c p 
(green dashed line) for the non interacting case, and c p (blue 
dot-dashed line) in the interacting case. Inset: behavior of 
the pressure as a function of the temperature, with and with- 
out interactions. Black dot-dashed line: pressure for the non 
interacting case. Red solid line: pressure for the interacting 
case. In all cases the density is fixed at n = 10 22 cm" 3 , and 
only first and second ionizations are included. 



tky anomaly discussed above on the temperature achiev- 
able in a process of adiabatic compression. We focus on 
this particular thermodynamic transformation because it 
plays a fundamental role in different atomic phenomena 
of great recent interest, such as sonoluminescence. In- 
deed, in sonoluminescence the (mostly noble) gas inside 
a bubble follows an isothermal transformation during the 
expansion of the bubble and an adiabatic compression 
during the bubble's collapse, the latter being the charac- 
terizing phase for the emission of light [l3j . In principle, 
kinetic Schottky anomalies in correspondence of ioniza- 
tion processes should have a sizeable effect on the col- 
lective heating and on the temperature achievable in an 
adiabatic compression. On the other hand, these types of 
transformations for gases in intermediate regimes of den- 
sity and temperature are typically described by the equa- 
tions of state of ideal gases supplemented by their van dcr 
Waals refinements to account for corrections due to in- 
teractions. In cither case, ionization processes cannot be 
handled properly, and Schottky-like behaviors cannot be 
predicted. We will thus compare the description of adia- 
batic compression obtained using the ideal-gas equation 
of state, the van der Waals one, and our self-consistent 
set based on the full microscopic description of atoms and 
ions. To set the framework, in the following we consider 
a quasi-static adiabatic compression of a sphere contain- 
ing one mole of atoms which is initially in thermal equi- 
librium under standard conditions at room temperature: 
P = latm, T = 300°A'. An adiabatic compression is 
performed until the density in the sphere becomes com- 
parable with the one present in a sonoluminescent bubble 
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FIG. 4: Main plot: behavior of the temperature inside a bub- 
ble containing a mole of atomic g eLS ELS Sb function of the radius 
of the bubble (expressed in unit of the ambient radius R a ) af- 
ter an adiabatic compression. The prediction of the ideal-gas 
equation of state (black solid line) is compared with those 
obtained from our set of state equations for non interacting 
Helium (blue dotted line), non interacting Argon (green dot- 
dashed line) and non interacting Xenon (red dashed line). In- 
set: total number of free electrons, as a function of the bubble 
radius, in the case of an adiabatic compression for the same 
atomic species considered in the main plot. 

during the last stage of its collapse. We investigate the 
predictions on the behavior of the temperature inside the 
bubble after the compression, according to the standard 
state equations (ideal or van der Waals) and to our set of 
state equations combining the description of the atomic 
interactions and the ionization processes with Lennard- 
Jones type potentials in a unified framework. 

The results of our comparison are reported in Fig. 0] 
Notice that when the radius of the sphere R is suffi- 
ciently large, the solutions obtained modeling the atomic 
mole as an ideal gas are indistinguishable from the mi- 
croscopic model Eqs.([T])-([5]), as it should be. However, 
when R becomes comparable to or smaller than 0.2 R a , 
where R a ~ 18.04 cm denotes the radius of the sphere at 
standard room temperature conditions, the central tem- 
peratures predicted by state equations Eqs.([T])-([5) are 
considerably lower than the ones obtained for the ideal 
gas. The lowering of the temperature is due to the fact 
that a part of the energy provided to the system by the 
compression is used up by some of the atoms into ion- 
ization processes. Xenon is the gas in which this effect 
is more relevant because its ionization energies are lower 
than the ones of Argon and Helium. Even if the num- 
ber of free electrons is low compared to the number of 
non-ionized atoms, and hence the system is only at the 
onset of the kinetic Schottky anomaly, the effect is so 
strong to halve the temperature as, e.g., in the Helium 
case at R = O.Q5R a . The correction is even stronger if 
we compare our results with the temperatures predicted 
using the van der Waals state equation. Indeed, the lat- 
ter induces an unphysical divergence in the temperature 
when the volume of the bubble approaches the co- volume 
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of the atomic species involved. Finally, we observe that 
for densities comparable with the ones typical in sono- 
luminescent bubbles, the electronic density is in good 
agreement with the most accurate experimental results 
on sonoluminescence [i~8j , according to which the plasma 
electron density in the bubble exceeds 10 21 cm -3 . 

Conclusion - We have introduced a set of state equa- 
tions for warm ionized gases, characterized by the fact 
that the interatomic, atom-ion, and interionic interac- 
tions are modeled by accurate microscopic potentials of 
the Lennard- Jones type, and the populations in the se- 
quential ionization processes are not fixed arbitrarily a 
priori but are obtained self-consistently using the Saha 
equation. We have analyzed in detail the case of noble 
gases, because of their importance in many physical phe- 
nomena such as sonoluminescence. However, our work 
can be extended straightforwardly also to different kinds 
of atoms, molecules, and mixtures of them by allowing for 
chemically bounded elements. Using this new approach 
we have shown that the noble gas atoms at low densi- 
ties exhibit a series of anomalies (fast increase followed 
by a fast decrease) in their specific heats, that represent 
the kinetic analogous for continuous energy spectra of 
the Schottky anomalies observed in the specific heats of 
solid-state systems with discrete and finite energy spec- 
tra. Moreover, we have shown that pressure and density 
of the free electrons exhibit a step-wise behavior in cor- 
respondence of the different ionization processes. The 
anomalies in specific heats become less evident at high 
densities and disappear when all the atoms are fully ion- 
ized. We have also shown that interactions affect rele- 
vantly both the pressure and the specific heat at fixed 
pressure c p while, on the contrary, the specific heat at 
fixed volume c v and the electron density n_i remain es- 
sentially unaffected. 

This novel type of kinetic Schottky anomalies predicted 
for the specific heats of noble gas atoms can play a very 
relevant role in many phenomena of actual interest, for 
instance sonoluminescence, in which collapsing bubbles 
close to their minimum radius are almost completely 
filled with noble gas atoms, reach temperatures of few 
tens of thousands of Kelvin, and densities of the order of 



the ones typical of solids |14| . In such phenomena, 
during the last stage of the collapse phase, the tempera- 
ture and the radial speed of the bubble depend crucially 
on c p and c„, so that their increase in correspondence of 
a kinetic Schottky anomaly can influence the dynamics 
of the bubble quite strongly. In the present paper we 
have not directly considered the problem of the full dy- 
namics of sonoluminescence, which will be analyzed in 
detail elsewhere, but we have only applied our general 
findings on the thermodynamics of warm ionized gases 



to the simple problem of the adiabatic compression of a 
sphere that reaches a final density comparable with the 
one of sonoluminescence during the emission of light, and 
we have predicted values of the electronic density for the 
gas in the sphere that are in very good agreement with 
the currently available experimental results. 
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